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Abstract 

We explore models in which weakly interacting massive particle (WIMP) dark matter annihilation is 
directly responsible for baryogenesis, thereby connecting dark matter with baryogenesis. We call this 
process "WIMPy baryogenesis". The dark matter relic density in these models, as with conventional 
WIMP models, is obtained with only order one couplings and TeV-scale masses according to the WIMP 
miracle. Thus, WIMPy baryogenesis models naturally accommodate weak-scale dark matter. Furthermore, 
an extension of the WIMP miracle simultaneously explains the observed baryon asymmetry and the correct 
dark matter abundance. The models we present have the further feature that they create the baryon number 
asymmetry at the weak scale, thereby avoiding the problems in some models of baryogenesis associated with 
high reheat temperatures in supersymmetric theories. Some of these models yield observable consequences 
in ongoing and future experiments. 
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1 Introduction 

Generally, baryogenesis and the establishment of the dark matter number density are treated as independent 
processes. With the notable exception of one class of models [7], it has largely been overlooked that models 
with symmetric, weakly interacting massive particle (WIMP) dark matter can connect dark matter physics 
with baryogenesis. We present a new mechanism that creates such a link and is based on a simple premise: if 
WIMP annihilation satisfies the Sakharov conditions, a non-zero baryon number asymmetry can be generated 
from dark matter annihilation, and in some instances, can account for the entire observed baryon asymmetry. 
We call this process WIMPy baryogenesis. Our models are distinct from models of asymmetric dark matter, 
which propose that dark matter and baryons have their origins in a common asymmetry. In our models, 
the energy densities of baryons and dark matter are more loosely linked but can accommodate the observed 
dark-matter-to-baryon ratio. 

[HElEllllElE] 

We list below the Sakharov conditions and how they are satisfied in WIMPy baryogenesis: 

1. Baryon number violation: WIMP annihilations violate baryon or lepton number. A preserved U(l) 
symmetry is allowed if the baryon asymmetry is balanced by a negative asymmetry in a decoupled sector 
that restores the net global symmetry. We have such a U(l) symmetry in most models we present. 

2. CP violation: WIMP couplings to Standard Model fields violate CP. 

3. Departure from thermal equilibrium: The cooling of the universe provides the necessary departure from 
thermal equilibrium. Net dark matter annihilation begins around temperatures T < todm, resulting in 
a small deviation of the dark matter number density from its equilibrium value. The annihilations can 
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Figure 1: Schematic diagrams showing the evolution of the asymmetry created by dark matter annihilation, 
(left) Model where asymmetry created in exotic antibaryons is sequestered in a sterile sector through baryon- 
number-conserving decays, (right) Model where asymmetry created in exotic antibaryons is converted into a 
Standard Model baryon asymmetry through baryon- number-violating decays. 



generate a baryon asymmetry that depends on the amount of dark matter annihilation occurring during 
washout freeze-out, which is comparable to the dark matter density at that time. 

We present models that satisfy all three Sakharov conditions, and that simultaneously generate the observed 
baryon asymmetry and WIMP relic density. In particular, we find that there exist successful models of WIMPy 
baryogenesis with 0(1) couplings and CP phases, and weak-scale masses for all new fields. This is an extension 
of the WIMP miracle to also include baryogenesis, although we show that the size of the generated asymmetry 
is sensitive to the parameters in the theory and can vary by several orders of magnitude from the observed 
asymmetry. 

Although WIMP annihilation can generate a baryon asymmetry, there are other processes that have the 
potential to wash out the asymmetry, and their freeze-out is crucial to create the observed baryon asymmetry. 
In our models, the two leading sources of washout are inverse annihilations of baryons into dark matter 
and baryon-to-antibaryon processes. Washout scatterings must be suppressed to generate a sizeable baryon 
asymmetry because, as we show in Section [21 any asymmetry generated prior to washout freeze-ou10 is rapidly 
damped away. After washout processes freeze out, dark matter annihilations can efheiently create a baryon 
asymmetry, and the final asymmetry depends on how much dark matter remains when washout scatterings 
freeze out. Washout freeze-out must occur before that of WIMP annihilation, at which point dark matter 
annihilation is no longer efficient and no sizeable asymmetry can be created. Thus, we find our central result: 
if washout processes freeze out before WIMP freeze-out, then a large baryon asymmetry may accumulate, and 
its final value is proportional to the WIMP abundance at the time that washout becomes inefficient. 

The early freeze-out of washout processes can occur for kinematic reasons. Inverse annihilations will be 
Boltzmann-suppressed for T < ttidm because the thermal baryon fields are no longer energetic enough to 
annihilate back into dark matter. Baryon-antibaryon scatterings, however, can remain rapid at temperatures 
well below ttidm- The only way to suppress baryon-to-antibaryon washout is if all washout processes involve a 
heavy exotic baryon field in the initial state. We illustrate this scenario in Figure [U showing how dark matter 
annihilates to Standard Model baryons plus an exotic baryon, as well as the possible decays of the exotic 
baryon (either through baryon-preserving or baryon- violating interactions). If this exotic field has a mass 
> mDM, its abundance is Boltzmann-suppressed at T < mDM and suppresses the washout rate. Meanwhile, 
dark matter annihilations are not kinematically allowed if the heavy baryon field has mass > 2todm, so the 
mass condition ttidm ^ "^-exotic baryon ^ 2todm is essential to generate a large baryon asymmetry through 
WIMPy baryogenesis. 

Dark matter annihilations generate a positive baryon asymmetry stored in Standard Model quarks along 
with an equal negative asymmetry stored in the exotic baryon field. It is important that the decays of the 

^The time of washout freeze-out is defined as when the rate of washout processes falls below the Hubble expansion rate. This 
is analogous to the freeze-out of WIMP annihilation. 
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exotic baryon do not eliminate the Standard Model baryon asymmetry. In models of WIMPy baryogenesis with 
a preserved U(l) baryon symmetry, the exotic baryon-number-carrying field is charged under an additional 
discrete symmetry, while Standard Model fields are uncharged, preventing the exotic baryon from decaying 
into Standard Model baryons and destroying the asymmetry. The heavy baryon-number-carrying field decays 
instead into light gauge singlet fields that are charged under the discrete symmetry and decoupled from 
Standard Model fields at temperatures below the scale of WIMPy baryogenesis. We also present a model 
where the exotic baryon decays to Standard Model quarks through baryon-numbcr-violating couplings, and 
such models manifestly satisfy the first Sakharov condition. 

For many years, the WIMP miracle - the fact that dark matter fields with weak-scale masses and annihila- 
tion cross sections give the correct dark matter thermal relic density - has been a compelling paradigm for dark 
matter model-building. WIMPy baryogenesis preserves the WIMP miracle while also offering an explanation 
for the observed baryon asymmetry. While WIMPy baryogenesis models do not predict the precise relation- 
ship between the dark matter and baryon number densities, natural models do restrict the baryon asymmetry 
to a range of about seven orders of magnitude (see Section [2]) , and the observed asymmetry is within this 
range. Since baryogenesis arises from WIMP annihilations, WIMPy baryogenesis is also necessarily connected 
to weak-scale physics. While we do not discuss an embedding of WIMPy baryogenesis in a particular solution 
of the hierarchy problem, we assume that whatever new physics lies at the weak scale stabilizes any scalar 
potentials in our theory and gives a natural explanation for their weak-scale masses. A consequence of this 
is that, with weak-scale masses, some of the new fields necessary for baryogenesis may give signals at future 
experiments. Additionally, the present-day dark matter is symmetric, leading to the possibility of the indirect 
detection of dark matter annihilations as in conventional WIMP models. This is in contrast with generic 
asymmetric dark matter models, in which the majority of dark matter annihilations ceased long before the 
present day, although it is also noteworthy that there do exist scenarios in which the symmetric component 
of dark matter is regenerated at late times, giving indirect detection signals for some asymmetric dark matter 
models [6]. 

A further advantage of this scenario is that bounds on the reheat temperature in supersymmetric models 
do not constrain WIMPy baryogenesis. Typical reheat temperature constraints come from overproduction of 
gravitinos and are in the range Trh S 10^ ~ 10^ GeV [5], and Trh is consequently below the scale required for 
conventional leptogenesis through the decay of heavy, Majorana right-handed neutrinos. Although low-scale 
mechanisms for baryogenesis are known, such as electroweak baryogenesis, WIMPy baryogenesis is a new way 
of generating the baryon asymmetry at T ^ TeV while satisfying the reheat bound. 

We discuss the general conditions for successful WIMPy baryogenesis in Section [21 finding that interactions 
washing out the baryon asymmetry must become ineffective prior to WIMP freeze-out in order to generate the 
observed asymmetry. In Section [31 we focus on a particular model where dark matter annihilates through a 
lepton-number-violating interaction and the asymmetry is subsequently transferred to baryons by sphalerons. 
Because sphaleron processes are only rapid in the unbroken electroweak phase, such baryogenesis must occur 
before the electroweak phase transition. We compute the dark matter relic density and baryon asymmetry, 
and find the range of masses and couplings that agrees with the observed densities of both. We also consider 
the implications of additional lepton-number-conserving dark matter annihilation channels. In Section [31 we 
consider models with WIMPs annihilating directly to quarks, where baryogenesis can occur over a wider range 
of temperatures because sphalerons are no longer needed to establish the baryon asymmetry. We discuss 
experimental constraints and possible signals for models of WIMPy baryogenesis in Section [5l Finally, we 
summarize in Section [6l 

2 General Analysis of WIMPy Baryogenesis 

To begin our discussion of WIMPy baryogenesis, we highlight in Section [2T] some of its general features and 
use an analytic approximation to determine the regimes in which baryogenesis is successful. Our central result 
is that the final baryon asymmetry from WIMPy baryogenesis is proportional to the dark matter density at 
the time when washout processes freeze out. This means that washout scatterings must freeze out at a time 
when the dark matter density was larger than or comparable to the observed final baryon asymmetry, and 
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that washout freeze-out occurs at such a time when one of the baryon-number-carrying products of WIMP 
annihilation is heavier than the dark matter mass, as described in the introduction. In Section [2.21 we estimate 
the magnitude of the baryon asymmetry for input parameters consistent with the WIMP miracle and show 
that it can lie within a range of approximately seven orders of magnitude. 



2.1 Boltzmann Equations and Solutions 

Wc consider a theory with dark matter species X whose annihilation violates baryon number, creating one 
quark (or anti-quark) along with other field(s) ipt. In this section we will not specify the precise interactions 
mediating dark matter annihilation in order to avoid specific model dependence. X can be either Majorana 
or Dirac, and the results derived below remain the same up to 0(1) multiplicative factors. All manifestations 
of the Boltzmann equation for baryon number evolution in WIMPy baryogenesis models have two important 
terms: one that describes the annihilation of dark matter and the consequent generation of a baryon asymmetry, 
and another that drives the asymmetry towards its equilibrium value of zero through baryon-number-violating 
washout scatterings. The full Boltzmann equations describing the evolution of the various particle abundances 
are model-dependent and can have many terms, which wc give explicitly for concrete models in Sections [3] and 
m However, in the models of interest to us, namely models where the asymmetry arises predominantly from 
WIMP annihilations, the overall dynamics arc well-described by the inclusion of only these terms. 

Consider the limit where WIMP annihilations are the dominant source of the baryon asymmetry and for 
which the asymmetry is small as observed. We derive the Boltzmann equations in terms of dimensionlcss 
quantities: the number density per comoving volume of field i, Yi = rii/s (s is the entropy density), and the 
temperature, which we express as a; = mx /T. The dark matter number density is denoted Yx and the baryon 
asymmetry is denoted Yas • The Yx evolution equation has one term that is important in all models of WIMPy 
baryogenesis, namely the conventional WIMP annihilation term that is proportional to the annihilation cross 
section (Tann and drives Yx to its equilibrium value. This term arises from both XX — > baryon processes and 
the inverse processes, baryons — > XX . The Yx Boltzmann equation with this term is 

^ = (-an„«) - [Y^r] , (1) 

ax xH[x) A J 

where H{x) is the Hubble scale. 

We neglect a back-reaction term in the Yx Boltzmann equation e s{x) {(Ta.auv)Y/\B{Y^Y / {2Y^ x H{x)), 
where e is the net baryon number created per dark matter annihilation and is a measure of the magnitude of 
CP- violation (it is defined more precisely in Section [3]). This term accounts for the modification of the inverse 
scattering rate of baryons into X when there is a baryon asymmetry. This approximation is valid because this 
term is small when Yab ^ 1, as is true in our universe (Yab ^ lO^^''). This simplification also decouples 
the equations for Yx and Yab, which makes it easier to get an approximate analytic solution for Yab- The 
equation ([1]) in this limit is the same as the familiar Boltzmann equation for conventional WIMPs. Yx is then 
obtained from the standard WIMP relic density calculation and is approximately inversely proportional to 
the annihilation cross section. 

The Boltzmann equation for the evolution of the baryon asymmetry has two important terms. In the 
first term, a baryon asymmetry is generated through X annihilations, and is proportional to e/2 x dYx/dx, 
which is the annihilation rate multiplied by the fractional asymmetry generated per annihilation. The factor of 
1/2 arises because the annihilation term in ([T]) includes the sum of annihilation into baryons and antibaryons, 
whereas the term generating the asymmetry includes the difference. The second term in the baryon asymmetry 
Boltzmann equation reduces the existing baryon asymmetry and is the washout term. It is proportional to Yab 
multiplied by the cross section of processes that eliminate the baryon asymmetry (Twashout- The Boltzmann 
equation is 

- "^"^ Va„„«) [y| - (y-)^] - ^ (.wa.hout.) ^ n ^r- (2) 



dx xH{x) ' L ' J xB{x) 2Y 
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The factor of Yab/2Y^ comes from the the fact that the chemical potential /iab for the baryon asymmetry 
can be written as fJ-As/T ~ Yab/'2Y^ [5]. We assume that all species except for X arc in equilibrium. There 



4 



are other terms that we have not included, such as washout terms proportional to i^oxotic-s ^ that come from 
scattering of baryon- number-carrying fields off dark matter fields. Typically, the suppression coming from the 
small value of Yx for a; S> 1 makes this term subdominant to other washout terms. In the models of Section 
[3]and[4l the i^oxotic-s term can be ignored without substantially affecting the numerical results. 

As expected, the Boltzmann equations show that the total baryon number is zero when all fields are in 
equilibrium and with the initial condition Yab = 0. A solution for Yab can be written in integral form in 
terms of the X density: 



YabIx) 



dx' 



esjx') 
x' H{x') 



^ , dYxjx') 

dx — — exp 

dx' 



[Yx^f] (^') cxp 



" dx" s{x") 



dx" s{x") 



2Y^ H{x") 



X'^washout'^, 



2Y^H{x") 



YlYrix") 



(o-washoutw) Yl_ Yi'^ix") 



(3) 



(4) 



Equation Q explicitly shows that YAsix) can be expressed in terms of a source term from dark matter anni- 
hilations, and an exponential term that attempts to erase any asymmetry generated by WIMP annihilations. 
The source term can be written as dYx/dx, as in [TU]. At T > mx, or a; < 1, WIMP annihilations are balanced 
by inverse scattering processes and dYx/dx ~ 0, meaning no asymmetry is generated according to (|4]). At 
^ Ij the expansion and cooling of the universe result in net WIMP annihilations (dYx/dx ^ 0), providing 
the departure from equilibrium necessary for baryogenesis. The net asymmetry at any x is sensitive to the 
rate of washout processes during the epoch of WIMP annihilations. 

The integrand in the exponent is the washout rate rwashout(a^) normalized to the Hubble scale H{x), 

Twashout (*^) ^{x^ 



H{x) 



2Y^ H{x) 



shout 



(5) 



Washout freezes out when FwashoutZ-ff < 1- In the limit where Fwashout is a rapidly decreasing function of x, 
^ takes a particularly simple form. This is true if, for example, the washout rate freezes out because rrii/T 
becomes large and yields an exponential suppression of Y^'^. In this case, we can model the exponential in (j4]) 
as a step function and obtain 

,dYx{x' 



Yab{^) « / 



dx' 



dx 



shout 



(6) 



where Xwashout = wx/TWashout is the point at which washout processes freeze out, and Yx(oo) is the late-time 
dark matter relic density. 

Equation ([6]) has a very clear physical interpretation: after washout scatterings freeze out, all subsequent 
WIMP annihilations generate a baryon asymmetry with efficiency e. This is why, according to Yab 
is proportional to e times the total number of WIMP annihilations that happen after Xwashout, which is 
Yx{x^a.s\io\it) — Yx{oo). The observed baryon asymmetry is Yab ~ 9 x 10~^^ [TT]. Since dark matter at late 
times satisfies the relation 

(5 GeV)rAs(oo) 



rx(oo) 



mx 



(7) 



we require that Yx{oo) < YasIoo) for weak-scale dark matter. Along with the requirement e < 1, ([7]) and ([6]) 
imply that Yx (^^washout ) 3> Yx(oo). In other words, the washout interactions must become ineffective prior 
to XX annihilation freeze-out in order to generate a sufficiently large baryon asymmetry through WIMPy 
baryogenesis. As an example of the numerical scales in WIMPy baryogenesis: for a WIMP of mass 1 TcV, 
Yx{oo) « 4 X 10"" and WIMP freeze-out happens at Xf,o. ~ 27, or T w 37 GcV. For e = 0.1, washout 
scatterings must freeze out at x^ashout ~ 20 or T ss 50 GeV. The final baryon asymmetry is proportional to 
the WIMP density at the time when washout ceases to be important, with yAB(oo) w 9 x 10"^"'^. 

For what parameters do we expect washout processes to freeze out prior to WIMP annihilation freeze-out? 
We compare Fwashout in (O to the corresponding rate of WIMP annihilation, which is [9] 

FwiMp(a:^) 2s{x) 



H{x) 



H{x) 



Yxix). 



(8) 



5 



We then find that 



r washout jx) ^ (gwashout^) Hi 



(9) 



rwiMp(a;) A{ai,nnv)Yp{x)Y^ 



This ratio must be less than one at the time of washout freeze-out for washout processes to freeze out prior 
to WIMP freeze-out. This can be realized if either of the following is true for every process washing out the 
baryon asymmetry: 



2. The baryon-nuniber-violating coupling is small so (cTwashout v) <C (cTann v). 

The second scenario is challenging to realize, because the same baryon-number-violating couplings appear 
in both the washout and annihilation cross sections, and (cTann v) is fixed by the dark matter relic density. 
Furthermore, as we show in Section [3.1.3[ suppressing the washout cross section also suppresses the fractional 
asymmetry generated per annihilation, e, and the resulting baryon asymmetry is typically too small. Therefore, 
we expect that viable models of WIMPy baryogenesis have at least one baryon-number-carrying field with mass 

2.2 Estimates of Baryon Asymmetry 

In this section, we derive an estimate of the baryon asymmetry generated by a WIMP dark matter candidate 
with mass mx ~ TeV, and we determine the size of the baryon energy density compared to the WIMP relic 
density. In the following, we assume for simplicity that a dark matter field X annihilates into a Standard 
Model quark Q plus an exotic baryon field (see Sections |3] and |4] for specific model details). We find that 
the baryon asymmetry depends strongly on the mass and is constrained to lie within a seven or eight 
order-of-magnitude window, with the observed baryon asymmetry within an order of magnitude of the upper 
limit. Therefore, WIMPy baryogenesis does not predict the value of the dark matter-baryon ratio, but neither 
is the relationship between the two energy densities completely arbitrary. 

To determine the range of baryon asymmetries obtained from WIMPy baryogenesis, we use the result from 
the last section that the final baryon asymmetry is proportional to the number of dark matter annihilations 
that occur after washout freeze-out, as shown in The largest possible asymmetry is generated when the 
exotic baryon field is heavy relative to dark matter (m^ > mx) so that washout processes freeze out while 
there is still a large dark matter abundance. To determine the upper bound on the asymmetry, we use the fact 
that < 2mx for WIMP annihilation to be allowed kinematically, and this limits how many dark matter 
particles can remain when washout freezes out. By contrast, the baryon asymmetry is small when washout 
processes turn off at a late time (m^ <C mx) after dark matter annihilation has frozen out. To calculate the 
lower bound, we determine the rate of residual dark matter annihilation after dark matter freezes out and use 
this to determine the size of the asymmetry. 

For both the upper and lower limits, we first calculate the allowed baryon asymmetry and then determine 
the corresponding dark matter-baryon ratio. In both scenarios, the baryon asymmetry depends on two time 
scales: the point of washout freeze-out, Xwashout = "T-x/^washout, and the point at which WIMP annihilation 
freezes out, a;ann = "ix/Tknn- 

Estimate of upper limit: We first estimate the upper limit of the baryon asymmetry generated within 
our framework, which occurs when m^ is heavy to suppress washout and is therefore also at the TeV-scale. 
Kinematically, dark matter annihilation occurs only if < 2mx, which bounds how early Xwashout can be 
relative to Xann- For a TeV-scale dark matter field, WIMP annihilation freezes out when the temperature is 
about 1/30 of its mass. Therefore Xwashout ~ Xannimx /m^) > 15. We also know that, when washout freezes 
out while WIMP annihilation is still active, (2;washout ) ^ ^(oo). We then obtain from ([6]): 



According to (jlOp , the asymmetry is independent of mx and depends only on the ratio x^ashout ~ m^ /mx , 
with a large asymmetry when m^ is comparable to or larger than mx ■ 



1 . One of the baryon states is heavier than dark matter so 




< 1. 



>AS(00) « - Ijf (Xwashout) < 2 Yxi^^) ~ £ X 10"'- 



(10) 
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To compare the baryon density to the dark matter energy density, recaU that the WIMP density changes 
httle after annihilation freezes out, and so Yx(oo) « ix(a^ann) ~ ix(30). We then find that 

_ mp,oton>AB(oo) £ Ix (a: washout ) / GcV \ 6 ^^^(15) GcV \ / £ A /^TcV\ 

Therefore, for a model with weak-scale mx, 0{1) couplings (in accordance with the WIMP miracle), and the 
loop-suppressed e ^ 10~^ as in we find that the energy density of baryons can be at most an order of 
magnitude larger than the energy density of dark matter. 

Estimate of lower limit: In deriving the upper bound, we assumed that saturated the bound < 2mx 
and we found that dark matter annihilation could generate the observed baryon asymmetry. When to^, ^ 2tox : 
washout processes remain in equilibrium until after dark matter freeze-out, and the asymmetry from WIMPy 
baryogenesis is too small to account for the observed asymmetry. We now estimate the full range of baryon 
asymmetries achieved in our models when ^ 2mx- In this case, the equilibrium number density of X is 
much smaller than the actual, frozen-out X abundance. As a result of this overabundance of X relative to its 
equilibrium value, some residual dark matter annihilations continue at late times, even though the annihilation 
rate is insufficient to appreciably change Yx after Xann- Such annihilations can, however, generate a small 
baryon asymmetry. According to ([5]), this asymmetry can be estimated by calculating Yx{oo) — yx(a;washout), 

where X washout ^ ^ann- 

To determine the asymmetry, we solve the Boltzmann equation ([l}, neglecting the subdominant term 
(y^'^)^ in equation ([T|). Furthermore, if XX annihilation is s-wave, then (iTann is approximately constant in 
the domain Xann < a; < ^washout- The only x-dcpendcnce comes from the factor 

■^(a^) '^(■^ann) '^a.nii (12) 



Integrating ([1]) from x = Xwashout to .t = oo gives 

-1^ / \ -ir / \ 2s(Xann)a^ann (f^ann '^^) (a^ann) 1 /-i o\ 

yx(a;washout) - Yx{od) « — r . (13) 

^ V'^ann ) -^washout 

Using the definition of Fwimp in ([HI), together with the fact that rwiMp(^ann) — H{xann)^ gives the simple 
result ^ ^ 

(-^washout 

)-rx(oo)« ^^Yx{x^nn)~ ^^Yx(oo). (14) 

•^washout -^washout 

Notice that for Xwashout > a;ann, Yx is constant at leading order from Xwashout to oo. Also as mentioned earlier, 
by assuming both X and have weak-scale masses and interactions, Xann/a^washout ^ m^,/mx- We can then 
obtain an estimate for the baryon asymmetry: 

Fab « Yxioo) « e ( ^) Yx{oo). (15) 

a^washout V ''^X / 

We see that the baryon asymmetry decreases linearly with m.0 when <C mx- 
The ratio of the baryon energy density to the dark matter energy density is 

_ 1 n-3 , "^'A \ f 



10-^X6^^ — . (16) 
ilx \mx J \ mx J 

If there is no large hierarchy in rax and (i.e. m^/mx ^ 0.1), and using our earlier estimate of e ~ 10~^ 
for 0(1) couplings that give the correct WIMP relic density, we find that ils/^x ^ 10"^. We emphasize, 
however, that even smaller asymmetries are possible if the imaginary parts of the couplings are tuned to be 
small or if there exist hierarchies in the masses of the new fields. 
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Figure 2: The evolution of the number density per comoving volume for field i {Yi) as a function of a; = mx /T . 
The numerical solutions shown here are based on the WIMPy leptogenesis model discussed in Section |3l where 
the dominant annihilation process is XX — Lif^ and the dominant washout is Lifj — > L^ip'^ . The input 
parameters are yx = 2.7, \l = 0.8, e = 0.2, mx = 3 TeV, and ms = 5 TeV. = 4 TeV gives the behavior 
when washout freezes out well before WIMP annihilation freezes out ("weak washout"), = 2 TeV gives 
the behavior when washout becomes ineffective subsequent to WIMP freeze-out ("strong washout"). 



Considering equations ([TT|) and ([TS]). we find that the expected range for the baryon-to-dark matter ratio 
in WIMPy baryogenesis is 

10-6 < % < 10, (17) 
ilx 

assuming 0(1) couplings, and weak-scale masses for all new fields, i.e. mx,rntj; ^ 0(0.1 — ITeV). The ob- 
served value of Hb/^x ~ 0.2 falls within this range, and thus WIMPy baryogenesis can account for the entire 
observed baryon asymmetry, but it does fall toward the upper end of the allowed region. 

To summarize, models of WIMPy baryogenesis predict a dark matter relic density inversely proportional to 
the WIMP annihilation cross section, as in conventional WIMP models, and a baryon asymmetry proportional 
to the dark matter density at the time when washout processes freeze out. In Figure [2l we illustrate the 
evolution of the dark matter abundance and the baryon asymmetry in one model of WIMPy baryogenesis for 
the two limiting washout cases. 



3 WIMP Annihilation to Leptons 
3.1 Model Overview 

We have discussed baryogenesis in the generalized sense of either the direct production of a baryon asymmetry 
through WIMP annihilation or leptogenesis, in which a lepton asymmetry is produced by WIMP annihilation 
and converted to a baryon asymmetry through sphalerons. In this section, we present a model of leptogenesis, 
where the lepton asymmetry is generated above the electroweak phase transition while sphalerons are still 
active. We discuss the field content and symmetries of the model Section [3.1.1[ and we calculate the efficiency 
of generating a lepton number asymmetry in Section 13.1.21 As we showed in Section [21 the final baryon 
asymmetry is determined by the time at which washout processes freeze out. We address washout in Section 
13.1.31 discussing the implications for the WIMPy leptogenesis parameter space. Finally, we give the Boltzmann 
equations in Section r3. 1.41 
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Table 1: Z4 charges of fields in the WIMPy leptogenesis model described by ([T^ . 



3.1.1 Field Content and Lagrangian 

Wc consider a simple model with the minimal ingredients for WIMPy leptogenesis. Dark matter consists of 
a pair of gauge singlet Dirac fcrmionsH X and X that annihilate to the Standard Model lepton doublet Li 
and new weak-scale fields ipi. X annihilates through weak-scale gauge singlet pseudoscalarto Sa- By gauge 
invariance, ipi has charge (2, 1/2) under the SU(2)l x U(l)y gauge interactions. The Lagrangian is 

C = £kin + 'Cmass - ^ (Axa^^ + >^'xa^^) + i >^L ai SaUlpi + h.C. (18) 

To satisfy the Sakharov conditions, dark matter annihilation must also violate CP, and Xlm must be complex. 
To have physical CP violation, there must be more than one scalar Sa so there is a relative phase in their 
amplitudes. XX annihilation can then generate an asymmetry in Li, and the lepton asymmetry is subsequently 
converted to a baryon asymmetry by sphalerons. Because the symmetry preserved by sphalerons is B — L, a 
negative lepton asymmetry must be generated to account for the observed positive baryon asymmetry. 

A positive lepton number asymmetry also accumulates in ipi , and it is important that this positive asymme- 
try does not erase the negative asymmetry in Standard Model leptons. In our model, ipi decays into light gauge 
singlets that are decoupled from Standard Model fields at low temperatures. The asymmetry produced in 
ip is therefore sequestered in a sterile sector and the Standard Model asymmetry persists to the present timcjfl. 
A Z4 symmetry, with charges in Table [TJ forbids other operators that allow ip to decay directly into Standard 
Model leptons, thus preventing the erasure of the Standard Model lepton asymmetry. The symmetry also 
makes dark matter stable. 

In the simplest model, ipi decays to Ui + H through the interaction 

AC = X, H^n^i', + h.C. (19) 

We assume that ipi is vectorlike with a partner ipi in order to more readily satisfy electroweak precision 
constraints. Its mass is restricted by the LEP bound > 100 GeV (see Section [5T^ . 

After electroweak symmetry breaking, ipi mixes with the sterile neutrino rii, and we must ensure that 
the sterile neutrino satisfies overclosure constraints. Since ip is Dirac, we also include ip when diagonalizing 
the mass matrix and find that there remains a massless eigenstate even after the Higgs condenses. This is 
good, because light, weakly interacting thermal relics (such as sterile neutrinos) with masses > 0{eV) would 
overdose the universe, n could have a Majorana mass m„ < eV and still satisfy observational constraints, 
since the light eigenstate would have a mass < eV as well, but we take n to be massless in our model. 

The Lagrangian ((18]) is also invariant under a U(l)'^ lepton flavor symmetry that prohibits flavor-changing 
neutral currents but allows flavor-dependent couplings. Li, rii, and ipi have charges -|-1, -1-1, and —1, respec- 
tively, under the U(l)i factor of the flavor symmetry. We assume that the only source of flavor-breaking in 
the low-energy theory is through the neutrino mass matrices, and this effect is very small. 

Majorana dark matter field X does not work in this case because X must carry a complex charge for the model to generate 
a non-zero lepton asymmetry, as we show later in this section. 

•^We consider pseudoscalars instead of scalars because they do not have a velocity-suppressed XX annihilation cross section. 
A scalar S with couplings to X that are CP-violating with a large imaginary part would work as well. 

*If > mg, then tp can decay into S + and wipe out the lepton asymmetry. However, S then subsequently decays into 
either L + H + n'^ or Lt -|_ }{* -|_ n (when XX — >■ Lip is kinematically allowed, S — >■ XX is kinematically forbidden because 
2m, X > rn^ > mg), and the difference in the rates of these decays generates another lepton asymmetry. If the efficiency of 
asymmetry generation from 5 decays is comparable to that from XX annihilations, the asymmetry is comparable to the original 
asymmetry created from XX — > L'0. 
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Figure 3: Diagrams of tree and loop contributions to S decay. The difference between these rates and their 
conjugates generates a lepton asymmetry. 




Figure 4: Diagrams of tree and loop contributions to the XX annihilation cross section. The difference between 
these rates and their conjugates generates a lepton asymmetry. 



3.1.2 Asymmetry generation 

In this model, a lepton asymmetry can in principle be generated through two processes: the more conventional 
process of Sa decay into Liipi and their conjugates (or directly into Li + Ui + H if > ms), and XX 
annihilation into the same final states. We show these in Figures [3] and [U respectively, assuming that decay 
and annihilation occur predominantly through the lightest scalar 5*1. Existing work discusses the relevant 
processes for generating a lepton asymmetry through 2 — >■ 2 scattering [12], although the authors consider 
only high-scale models (T > 10^ GeV) with qualitatively different features than WIMPy Icptogenesis. CP- 
violating phases in our model appear in the interference between tree-level and one-loop diagrams. We define 
asymmetry factors for the decay of the lightest scalar 5*1 and for WIMP annihilations, respectively, in the 
manner of conventional leptogenesis: 



ajXX ^ yj,L,) + ajXX ^ V>,L,) - ajXX ^ v4l|) - ajXX ^ 



ti gives the fractional asymmetry generated per Si decay, while £2 gives the fractional asymmetry generated 
per XX annihilation. The precise values of ei, e2 in this case depend on the masses msa and the couplings 

To reduce the number of arbitrary parameters in our analysis, we make the following assumptions: 

• Dark matter annihilation occurs dominantly to only one flavor of lepton, and the couplings of all other 
leptons to Sa are zero. The non-zero couplings of the single lepton flavor are denoted \La- 



Dark matter annihilation and washout occur mostly through the lightest scalar, ^i, and we consider the 
rates of only these processes in our analysis. For concreteness, we require that the corresponding cross 
sections with intermediate ^2 to be less than 20% of the corresponding cross sections with ^i, giving 
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roughly 



^4 

™S2 



< 



™S2 



r; 4 ' 

5m|i 



(22) 
(23) 



We also assume that TO51 <^ ms2, so that the loop integrals in ei and £2 can be put in a simple analytic 
form ([M)) . 

• The physical CP phases are large i.e. Im(a) w a where a is some product of couplings appearing in 
scattering and decay amplitudes. 

None of these assumptions are required by WIMPy leptogenesis, and we make them only to simplify the analysis 
and its interpretation. Relaxing these assumptions would introduce much complexity into the Boltzmann 
equations while giving qualitatively similar results. The phenomenology does, however, depend to some extent 
on the flavor of leptons to which dark matter predominantly annihilates (see Section [5] for details). With the 
above assumption^. 



1_ lm{Xl,Xl%) {2mxf 



'52 



\2mxJ \2mxJ \2mx J 



(24) 



The expression for ei is the same but with 2mx — > TO51. Since we are most interested in the asymmetry 
from annihilation, £2 is the relevant parameter for WIMPy baryogenesis and we denote its asymmetry factor 
by e = e2. e is suppressed by 1/771^2 from the S2 propagator, and is proportional to (2mx)^ because the 
momentum flowing through the Si propagator in XX annihilation is -s/l = 2mx +0{T)^ where T mx, msi 
at freeze-out. Note that ()24|) vanishes when = 2tox, at which point the particles in the loop cannot go on 
shell and there is no imaginary part of the amplitude (and, hence, no CP violation) . 

Using ((22|) . (p4)) . and the assumption of large CP phases, we can bound e from above: 



< 



'SI 



\2mxJ \2mxJ \2mx J 



(25) 



We treat e as a free parameter, subject to (|25p . and we can now express all rates and cross sections in terms 
of Ax = Ajfi, Al = Ail, e, mx, m^, and ms = msi- 

We have assumed that the lepton asymmetry from XX annihilations dominates over that from S decays. 
We find that this assumption is true whenever mx < ms- Since the asymmetry is proportional to the number 
density of X or S at the time of washout freeze-out, the ratio of asymmetry from decay vs. annihilation is the 
same as the ratio of the number of S particles to the number of X particles at the time of washout freeze-out. 
The assumption of annihilation-dominated asymmetry is therefore equivalent to mx < ms- 



3.1.3 Washout 

As we demonstrated in Section [2l the final baryon asymmetry depends on the time of washout freeze-out. 
We now discuss the implications for WIMPy leptogenesis, finding that we need > mx for successful 
WIMPy leptogenesis. We show the lepton number washout processes in Figure [5l They include inverse 
annihilations, lepton — > antilepton scatterings, and ipX — ^ L^X processes. The dominant washout is typically 

^This expression is derived in the narrow-width approximation. For TeV WIMPs, Xx ^ 1 is often necessary to obtain the correct 
dark matter rehc abundance, which may lead to Fgi, ~ ^S- When 2mx > fT^si, •S'l is kinematically forbidden from decaying 
into XX and the Si width is narrow (because typically Aj^ J; !)■ may be broad, but the imaginary part of its self-energy 
correction Imn(p^) as substituted into Ij24|l must be evaluated at = 4m^ <C ^"32 ^^'^ satisfies Imn(4m^) ~ <g; m^^. 

Similarly, if 2mx < msi, then the partial width of Si to X can be very large for 2mx ^ '"si i but once again Imll is evaluated in 
l|24|l as Imn(4m^) 4m^ < m^j^. Therefore, the narrow- width approximation holds true to a degree sufficient for our purposes. 
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Figure 5: Diagrams leading to washout of the lepton number from (top row) s-channel and (bottom row) 
i-channel scatterings. 



from Lip — > L^ip^ scatterings, because the inverse annihilation Lij] — > XX is kinematically suppressed for 
T < mx and ipX — > L^X gets more Boltzmann suppression. Applying ([5]) for the specific model of WIMPy 
leptogenesis, the washout rate is proportional to 

rwashout(.T) « 4? (^L.^^Lt^t V) Yl^Y^'^ix), (26) 

where s(x) is the entropy density at x. Washout freezes out when its rate is about equal to the Hubble scale, 
rwashout(a;washout) ~ (a^washout ) ■ Twashout (aJwashout ) cau bc smaU for one of two reasons: 

1. fTi-i/) ^ "nix so that ^^'^(xwashout) is Boltzmann-supprcssed while dark matter is annihilating. 

2- i'^Lip^L^ip'' ^) is small relative to the annihilation cross section so that washout freezes out before anni- 
hilation. The washout cross section can be small if <C 1. 

One of these two conditions must hold for each washout process. We find that option ^1 leads to viable 
WIMPy leptogenesis. Option #2, on the other hand, does not give a large asymmetry. According to 
the asymmetry efficiency factor e is also suppressed when Al <C 1, and the potential gain in the baryon 
asymmetry from early washout freeze-out in option ^2 is offset because leptogenesis occurs less efficiently. 
Therefore, > mx is generally required to generate the observed baryon asymmetry. 

When is much larger than mx (we find this is typically true for m,/, > 2mx), the exponential suppression 
of is so large that 3 — > 3 scatterings oi LnH — > H* become important (see Figure [S]). This region 

is, however, kinematically inaccessible in WIMPy leptogenesis since 2mx > fn^ for efficient annihilation to 
occur, and we neglect 3 — >■ 3 processes. 

3.1.4 Boltzmann equations 

We consider the evolution of a single component La, where a is a gauge index (flavor indices are suppressed 
since we consider only one flavor). We deflne Bil (Brx) as the total branching fraction of S into leptons 
(X), with Br^ -f- Brx = 1. Also, ^ = 1 + fJ-^^/fJ^L^, where fj, are chemical potentials, and 77 is defined as 
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Figure 6: 3 — > 3 washout process that can dominate over 2 — !• 2 scattering when T <^ m^. 



the amount of La asymmetry generated for each La directly created/annihilated (this accounts for the fact 
that the asymmetry spreads among ah baryons and Icptons by rapid thermahzation) . We calculated the cross 
sections and widths analytically and checked them with CalcHEP [13]. 

The Boltzmann equations describing the evolution of the various particle species and the asymmetry in 
one of the components of the La doublet are 



Hjmx) dYx 
X dx 



^rx{rs)Yp ( ^ ) + Brx(rs> {Ys - Br^ Y^) - e 



^Yal 
2K 



-BrxBvL{rs)YP; 



H{mx) dYs 
dx 



^ = ~{Ts)Ys + {Ts)Yp 



Brz, + Brx 



Yx 

X 



(27) 
(28) 



H{mx) dY^L^ 



XT] 



dx 



1 - 2BrL - Brx 



1 + 



y2 



+ 2se(axx«i„V'a [Yx - {Y''^f\ 



+ 



iY^ 
4Kv 



S {(JxX^ha.ila. ")(ix'')^ + 2S[(0- 



, t «>yxyr + 2s (a , 



BTL{rs)Yp (e'Bri+Brx). 



We assume that all abundances are in thermal equilibrium at a; = 1 and that all fields remain in equilibrium 
except for S and X . In the evolution of the scalar S, we only include the decay terms, as they dominate over 
SS annihilation for T <^ ms- 

To determine the relationship between /i^/,^ and /iL„ , we assume that all abundances other than S, X and 
ALa are in thermal equilibrium and that all processes except those involving S are in chemical equilibrium. 
We also take sphaleron processes to be in equilibrium. The non-S* couplings in (jlSp distribute the L and ip 
asymmetries among the light fields. Solving the chemical potential relations gives 



16+12nXK'^^ 
3 + 12nXK'^^ ' 

2i7 + 28n;ynll) 
79 + 355 n^/ nil ' 



(30) 
(31) 



(29) 



The precise values of ^ and ry are model-dependent. 
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The assumption of thermal equihbrium for is consistent provided the decay width Fy, > H(Ticp), which 
in our model constrains 

^2 > i6^g(riep) ^ &0V9:t,% ^ ^^^^ 

For Ticp = 100 GcV and m^i ~ 2 TcV, this gives Ai > 3 x 10"^. This is not a very stringent requirement, 
since this value is smaller than any of the Standard Model Yukawa couplings. 

The final lepton asymmetry is also determined by the chemical potential relations. The relation between 
the total lepton asymmetry ALtot and the asymmetry in a single component of the doublet field ALq as 
determined by equation ([^5]) is 



'5^ALtot - ^ , oo,cq°/^oq^° >ALa- (33) 



The final baryon asymmetry Iab follows from the sphaleron chemical potential relations, and is 



YM^) ^ -4 ( 1 y-L... = -4y.,„. (34) 



In the limit a; — >■ oo, the ratio of total baryon to lepton number reduces to the same expression as conventional 
leptogenesis [Hj . 

The total dark matter relic abundance is 



yDM(oo) = yx(oo) + y^(oo) = 2yjf(oo). (35) 



3.2 Numerical Results 



There are six free parameters in our model: three masses {ms, mx, and m^) and three dimensionless param- 
eters (Ax, Al, and e). To determine over what range of parameters WIMPy leptogenesis can be successful, we 
perform scans over two parameters at a time while holding others fixed. In particular, we are interested to see 
what range of masses is allowed, and if any tuning of the mass and coupling constant relations is necessary to 
generate the correct baryon asymmetry and WIMP relic density. 

Range of allowed masses: We hold ms fixed and determine for which mx and masses there exists some 
perturbative couplings that give the observed dark matter density and baryon asymmetry. We place no other 
restrictions on the couplings. If m^, > ms, we assume that the S width is dominated by the three-body 
decay S — > L H n. We show in Figure [7] the masses that give rise to successful WIMPy leptogenesis. The 
viable ip masses satisfy m^ (1 — 2)mx, while there is no correlation between mx and ms as long as 
mx < ms- For smaller values of m^/mx, the Boltzmann suppression of the washout rate is insufficient to 
generate the observed baryon asymmetry, while > 2mx is not allowed because dark matter annihilation 
is not kinematically allowed and because the asymmetry efficiency e is zero {CP violation is zero if L and ip 
cannot go on-shell in the dark matter annihilation loop diagrams). 

The lower boundary of the allowed region has a meandering shape around m^ « ms- The reason is that s- 
channel washout processes have a resonant enhancement in this region, leading to a smaller baryon asymmetry 
and a restricted parameter space. Above resonance, t-channel washout processes are also important, explaining 
why the bend in the curve is centered at m^ slightly larger than ms- 

The principal reason that it is difficult to generate a large baryon asymmetry is because the efficiency of 
asymmetry generation e is tied to the washout cross section through its dependence on X\ in p5|) . A large 
asymmetry can only be generated when washout effects are also large, limiting how much of an asymmetry 
can be generated. The viable parameter space is larger if ([25]) can be relaxed, as is the case when 5*1 and 5*2 
are nearly degenerate and the asymmetry is resonantly enhanced, but this is not a required feature of WIMPy 
leptogenesis. 

In leptogenesis, the asymmetry must be generated prior to the electroweak phase transition, at which point 
sphalerons decouple and the conversion of a lepton asymmetry into a baryon asymmetry ceases. Since WIMPy 
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Figure 7: Regions in the mx-m^ plane with the correct WIMP rehc density and baryon asymmetry from 
WIMPy leptogenesis, with ms = 5 TeV and some choice of pcrturbative couphngs. The masses giving both 
observed abundances are shown in blue (middle stripe). We plot the ratios mx/rns, m^/ms to show the 
relationship between the X and masses and the mediator scale ms- The excluded regions are shown in red: 
the upper region is not viable because 2mx < and the thermal annihilation cross section is Boltzmann- 
suppressed, while the lower region has too large to prevent rapid washout of the asymmetry. The dashed 
line indicates the lower boundary of allowed mx and m^; below the line, the electroweak phase transition 
occurs before the baryon asymmetry is large enough to account for the observed value. For mx/rns > 1, the 
asymmetry is dominated by S decay. 



leptogenesis is a weak-scale model, the timing of the asymmetry generation relative to the phase transition 
is important. To illustrate this, we computed the critical temperature Tc of the phase transition assuming a 
Standard Model Higgs with mass m^ = 120 GeV, and we required that the baryon asymmetry at Tc be equal 
to the observed asymmetrjf^. This typically yields a much smaller baryon asymmetry than lepton asymmetry 
at late times because the baryon asymmetry stops accumulating at . Accounting for the effects of the phase 
transition, the allowed region is above the dashed line in Figure [71 If the phase transition is modified by 
additional Higgs fields or other new physics, then this boundary line changes. 

Range of allowed couplings: We choose representative values of the masses, with ms ~ 5 TeV for all cases, 
and mx and m^, chosen in the middle of the allowed bands in Figure [T] For one set of parameters, dark matter 
annihilates above the S resonance, with parameters mx = 4.25 TeV, m^ — 7.5 TeV, and |e| — 0.075, and we 
determine the dark matter relic abundance and baryon asymmetry as functions of the two couplings. We also 
study XX annihilation below resonance, with mx = 1-5 TeV, m^ = 2.25 TeV, and |e| = 0.0075. We plot the 
results in Figure [8] as contours of constant relic density and baryon asymmetry. We focus on the ratio Xl/^x 
because we are interested in seeing if any relationship between these two theoretically unrelated quantities is 
required to obtain a particular relic abundance and asymmetry. In both cases shown, WIMPy leptogenesis 
gives the correct dark matter relic abundance and asymmetry when both couplings are 0(1). Thus, a perfectly 
natural choice of couplings, and the very same couplings that satisfy the WIMP miracle, can also generate 
the correct baryon asymmetry if CP phases are large! Specifically, with mx ~ 4.25 TeV, m,^ = 7.5 TeV, 
and |e| ~ 0.075, the observational constraints are satisfied with Xx = 2.7 and Xl = 5.7; with mx = 1-5 TeV, 
m^ = 2.25 TeV, and |e| = 0.0075, the couplings are Ax = 2.8, Xl = 2.5. 

^Since the Standard Model phase transition is of second order, sphalerons do not suddenly shut off, and a more proper treatment 
would account for the gradual departure from equilibrium of the sphaleron effects. Since the asymmetry is typically generated 
over a very short time period (we find numerically that it is on the order of Ax ~ 2 — 3 or AT ~ 5 — 10 GeV), the dynamics of 
sphaleron shut-off are largely irrelevant and the most important factor is the rate of L — > B transfer at the washout freeze-out 

time 3^washout ■ 
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Figure 8: Dark matter relic density (solid lines) and baryon asymmetry (dotted lines) as functions of couplings 
Ax and Xl/^x- We consider two sets of massses, with ms = 5 TeV for both: (left) mx = 4.25 TeV, = 7.5 
TeV, and |e| = 0.075; (right) mx = 1.5 TcV, m^, = 2.25 TeV, and |e| = 0.0075. The asymmetry contours 
are, from top to bottom: (left) Yab = 5 x 10~^^, 8.5 x 10~^^ (observed asymmetry), and 3 x 10~^°; (right) 
Yab = 3 X 10"^^, 8.5 X 10~^^ (observed asymmetry), and 3 x 10^^". The dark matter abundances are printed 
on the plots. In the shaded regions, the numerical value of e is not consistent with our assumptions according 
to the bound ((25)) . 



In deriving our results, we assumed that dark matter only annihilates through lepton-number-violating 
interactions. In a more general model, dark matter may also have lepton-number-preserving interactions that 
contribute to the total annihilation cross section. We parameterize this possibility with the quantity 

^^ (ax^^anythingt^) (36) 

{crxx^L'4, V) 

When a > 1, the asymmetry generated by WIMPy leptogenesis is smaller, because only 1/a of dark matter 
annihilations proceed through lepton-number-violating couplings and can create an asymmetrj0. As a result, 
the viable parameter space for WIMPy leptogenesis is reduced. In Figure [U we show the masses mx , Wi/i 
giving successful WIMPy leptogenesis with a > 1. In particular, wc find that the lepton asymmetry from 
WIMPy leptogenesis is too small in regions with large washout (m.0 ^ ms, where wash-out scattering is on- 
shell). While WIMPy leptogenesis is possible with lepton-preserving annihilation channels, m^j; lies in a more 
restricted region when a > 1. 

To summarize, we have presented a model of WIMPy leptogenesis where the WIMP miracle has been 
extended to the WIMPy baryogenesis miracle: the correct baryon asymmetry and WIMP relic density can be 
generated simultaneously with TeV-scale masses and C(l) couplings. We find that, depending on the ratio 
m^/mx and the Yukawa couplings, larger and smaller asymmetries are also possible over a range of about 
seven orders of magnitude. Generating the observed baryon asymmetry does require some correlation between 
mx and m^ , which may be explained if the masses have some common dynamical origin. For mx lighter than 
about 1 TeV, sphalcrons decouple in the middle of asymmetry generation and the resulting baryon asymmetry 
is typically smaller than the observed Yab- 



'^When a > 1, the WIMP annihilation cross section is also larger, and Xx, arc smaller to give the same WIMP relic density. 
As discussed in Section 13.1.31 however, decreasing the couplings results in both a smaller washout rate and a smaller efficiency of 
generating an asymmetry. These two effects counteract one another, and the change in couplings for a > 1 does not substantially 
affect the asymmetry. 
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ms = 5 TeV m,; = 5 TeV 




mx I ms mx I ms 

Figure 9: Regions in the mx-m^, plane for viable WIMPy leptogenesis with additional lepton-number- 
preserving dark matter annihilation modes: (left) a = 2, (right) a = 3. The masses giving the correct 
dark matter density and baryon asymmetry for some choice of perturbative coupling are shown in blue (mid- 
dle stripe). As the lepton- number-preserving annihilation cross section increases, the efficiency of asymmetry 
generation drops, and the marginal regions of parameter space become inaccessible, particularly the enhanced 
washout region ~ ms- The descriptions of the regions on the plot are the same as those in Figure [71 



4 WIMP Annihilation to Quarks 
4.1 Model Overview 

If the final products of dark matter annihilation arc quarks, WIMP annihilation can directly generate a baryon 
number asymmetry. The lower bound of ^ TeV on mx in WIMPy leptogenesis (the dashed line in Figure [71) 
docs not apply when WIMPs annihilate to quarks, since the production of baryon number no longer depends 
on efficient sphaleron interactions. Just as the leptogenesis model included new weakly charged vectorlike 
doublets, this model requires new vectorlikc colored states to couple to quarks. Such states can be pair- 
produced at the LHC, leading to much stronger constraints and better detection prospects, which we discuss 
in Section [Ol 

The model content is similar to the leptogenesis model discussed in Section [3l vectorlike gauge singlet 
dark matter X and X, singlet pseudoscalars Sa, and vectorlike exotic quark color triplets ipi and ipi. The 
Lagrangian is 

C = £kin + Anass - ^ {>'XaX'^ + >^'xa^^) + i>^B a SaUlp. (37) 

A baryon asymmetry can be generated in u along with an equal negative baryon asymmetry in ip. ip must 
decay, because it would otherwise overdose the universe and violate bounds on stable colored particles. 

The negative baryon asymmetry in tp must not destroy the positive baryon asymmetry in u when it decays. 
This can happen in two ways: 

1. "0 decays into a sector decoupled from Standard Model quarks at low energies. Total baryon number is 
preserved, but the negative baryon number carried by ip is sequestered from quarks at late times and 
does not eliminate the Standard Model asymmetry. 

2. -0 decays into Standard Model quarks through baryon-number-violating couplings. The final baryon 
asymmetry is different from the asymmetry created initially in u from WIMP annihilations because ip 
decays give an additional contribution to the baryon asymmetry. 

We now implement each of the above scenarios. 

1. decays to light, baryon-number-carrying singlets rii plus Standard Model antiquarks. It can do so 
through a colored scalar (p with Standard Model gauge representation (3, 1, —1/3). The additional terms 
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Tabic 2: Z4 charges of fields in models with WIMP annihilation to quarks. 



in the Lagrangian are 

A/: = A, dz 4>* + A', 4) d, n, + h.c. (38) 

A Z4 symmetry prevents -0 from decaying directly into Standard Model quarks through a QHrjj term 
and eliminating the baryon asymmetry. We show the charges in Table [5) This Lagrangian has a U(l)'^ 
flavor symmetry and satisfies all quark flavor constraints. In particular, Uj. tpi, and have charges — 1, 
+1, and —3, respectively, under the 11(1); factor of the flavor symmetry, cj) has charge —2 under all U(l) 
flavor symmetries. We assume that the only sources of flavor violation are the Standard Model Yukawa 
matrices. 

2. In this scenario, tjji also decays to two antiquarks plus a singlet n, but n is a Majorana fcrmion that does 
not carry any charge. Baryon number is now explicitly violated, and dark matter annihilations generate 
— 1 unit of baryon number for each + u produced from dark matter annihilations (because ip — > ddn). 
There is a new colored scalar di in the (3, 1, —1/3) representation of the Standard Model gauge group 
that mediates ip decays. The additional terms in the Lagrangian are 

AC = X e'^*^ rpi dj 4 + A ■ di d, n + h.c. (39) 

A Z4 symmetry, with charges given in Table [21 prevents from decaying to Standard Model quarks 
through other interactions that would destroy the baryon asymmetry. This Lagrangian can be naturally 
realized in supersymmetric models, where d is the down squark and n is the neutralino, although this 
is not the only possible realization of this scenario. (|39p has a U(3) flavor symmetry, which is the 
diagonal subgroup of the full U(3)„ x U(3)rf flavor group. The quark, ipi and di flelds transform in 
the fundamental of U(3). The Yukawa couplings between Su4> in (j37p have a flavor-independent piece 
and a flavor-dependent piece proportional to the up Yukawa matrix Y^, consistent with minimal flavor 
violation. 

In both scenarios, ip decays to a singlet n plus quarks. Operators allowing ^p to decay entirely to quarks 
(such as (p*du or d*du) are forbidden by the Z4 symmetry. The symmetry also ensures the stability of 
dark matter and of the proton. The proton is stable provided m-p < 2mx,'ms because baryons have charge 
(— 1)'^ ~ —1 and can never decay into the lighter meson and lepton flelds, which are uncharged under the Z4. 

The Z4 symmetry in principle allows neutral baryons to oscillate into one another. For scenario #1, the 
generalized baryon number symmetry prohibits neutron-antineutron oscillation. In scenario ^2, the baryon- 
number- violating term is antisymmetric in flavor indices, and the dominant contribution to neutron-antineutron 
mixing involves loops of W bosons and off-diagonal CKM matrix element insertions Vhd and Vsd- Since the 
bound on the neutron-antineutron oscillation operator c/A^{udd)'^ is A > 10— 100 TeV for c = 1 [15], the loop- 
and CKM-suppression is sufficient to lower the oscillation rate well below current constraints for TO5, to^, toj ^ 
TeV and 0(1 ) couphngs. 

The Boltzmann equations for WIMP annihilation to quarks are changed only by group theory factors from 
the corresponding equations for leptons. Similarly, the chemical potentials relations arc modified to reflect the 
new interactions ([38|l or ([39]) . 

As with WIMPy leptogcnesis, we assume that tp is in equilibrium, and this places constraints on the 
couplings through which it decays. In scenario #1, we considered the decay of ip according to interactions given 
in equation ([55)) . If > m^, the tp decay is two-body and the constraint (|32p applies, giving Ai > 6 x 10~*. 
If my, < m^, Ip undergoes a three-body decay to ddn, and the constraint is 
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Figure 10: Regions in the mx-m^ plane with the correct WIMP rehc density and baryon asymmetry from 
WIMPy baryogenesis, with (left) ms = 1.5 TeV and (right) ms = 5 TeV, and any choice of pertm^bativc 
couplings. The masses giving both observed abundances are shown in blue (middle stripe). The descriptions 
of the regions on the plot are the same as those in Figure [7l 



For Tiep = 50 GcV and m,^ = 1 TeV, we find that the constraint on the geometric mean of the couplings is 

The constraints on scenario #2 are comparable. 
4.2 Numerical Results 

Because of its similarity to leptogenesis, we use the quark flavor structure of scenario #1 in our analysis, since it 
allows for a more direct comparison of numerical results in both cases. We consider two scenarios: 7715 = 5 TeV, 
to compare the results for quarks with that for Icptons, and ms = 1.5 TcV, because dark matter can be much 
lighter than in WIMPy leptogenesis since there are no constraints from sphalcron decoupling. For simplicity, 
we consider sphalcrons to be out of equilibrium for the 1.5 TcV case and in equilibrium for the 5 TeV case to 
avoid considering sphalcron decoupling effects, although the calculation can be easily extended to include them. 

Range of allowed masses: We show the range of allowed masses in Figure [TUl Gluino searches at the LHC 
constrain this scenario (see Section r5.2p . in contrast with the leptogenesis model, for which the entire parameter 
space is unconstrained by collider searches. This is particularly true for ms = 1.5 TeV, where LHC searches 
will cover almost the entire parameter space for dark matter annihilation to quarks during the 14 TeV run. 
The WIMP mass is already constrained to be mx ^ 295 GeV by the gluino bound discussed in Section [5?2] 
along with the kinematic requirement that 2mx > m^. This is true independent of all other parameters. 

In both cases, the results are qualitatively similar to leptogenesis. With WIMP annihilation to quarks, the 
annihilation and washout cross sections arc enhanced because the final states arc charged under SU(3)c. As 
a result, the baryon asymmetry is suppressed by the increased washout rate. This is partially offset by the 
fact that the self-energy contribution to e is enhanced by a group theory factor as well. With mx = 0.9 TeV, 
the parameter space is actually larger than for mx = 4.25 TeV or WIMPy leptogenesis. Since in this case, 
sphalerons no longer inter-convert baryon and lepton number, the asymmetry created in quarks is distributed 
among fewer fields, enhancing the asymmetry. 

Range of allowed couplings: We do the same analysis as in Section 13.21 To compare the results of direct 
baryon asymmetry production with those for WIMPy leptogenesis, we choose a set of parameters used in the 
leptogenesis analysis: mx ~ 4.25 TcV, m^ = 7.25 TeV, ms = 5 TeV, and |e| = 0.075. We also consider a 
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Ax Ax 

Figure 11: Dark matter relic density (solid lines) and baryon asymmetry (dotted lines) as functions of couplings 
\x and As /Ax- We consider two sets of massses: (left), mx ~ 4.25 TeV, = 7.25 TeV, ms = 5 TcV, 
and |e| = 0.075. The asymmetry contours are, from top to bottom: Yab = 4 x 10"^"'^, 8.5 x 10"^"'^ (observed 
asymmetry), and 1.5 x . (right) mx = 0.9 TcV, m^, 1.2 TcV, ms ^ 1.5 TeV, and |e| = 0.075; the 

asymmetry contours are, from top to bottom: Yab = 5 x 10^^^, 8.5 x 10~^^ (observed asymmetry), and 
3 X 10^^°. In the shaded regions, e is not consistent with our assumptions according to the bound 



corresponding point with light S and broken phase chemical potential relations: mx = 0-9 TeV, m^, = 1.2 TcV, 
ms = 1.5 TeV, and |e| = 0.075. We show the results in Figure ITT] The parameter points giving the correct 
dark matter density and baryon asymmetry are Ax = 2.7, \b = 4.5 for mx = 4.25 TeV, and Ax = 0.22, 
As 2.8 for mx = 0.9 TeV. 

5 Experimental Constraints and Detection Prospects 

In this section, we survey the possible experimental constraints and signals for models of WIMPy baryogencsis, 
considering both annihilation to leptons and annihilation to quarks. For WIMP annihilation to Icptons, the 
experimental bounds on mx and m^ are too weak to constrain leptogenesis because mx , m^ > TeV are 
required to generate a sufficiently large baryon asymmetry. The prospects are better for WIMP annihilation 
to quarks, which predicts signals at indirect and direct detection experiments, as well as at the LHC. We first 
give a preview of our results in Table [3l 





Annihilation to leptons 


Annihilation to quarks 


Direct detection 




mx ^ 5 TcV for some parameters * 

(fTX-nuclcon ^ lO^^^ " 10-^4 cm^) 


Indirect detection 
(antideuterons) 


mx < 200 GeV 


mx < 1 TeV 


Colliders 


m^ < few hundred GeV, 
possible improvements 
with targeted searches 


m.0 < 1.44 TeV with 
100 fb~^ LHC (14 TeV) 


EDM 







Table 3: Search reach for minimal models of WIMPy baryogenesis/leptogenesis in current and near- future 
experiments. ' ' indicates no signal in that search channel. 
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Figure 12: Feynnian diagrams for dark matter scattering off (left) nucleons and (right) electrons in direct 
detection experiments for WIMPy leptogenesis. 




Figure 13: Feynman diagrams for dark matter scattering off nucleons in direct detection experiments when 
WIMPs annihilate to quarks: (left) standard signal and (right) inelastic induced nucleon decay. 



We provide details for each class of experiment in the following sections. 

5.1 Dark Matter Detection 
5.1.1 Direct Detection 

Dark matter direct detection experiments are typically important probes of weak-scale dark matter models. 
As we show in this section, however, only WIMPy baryogenesis with dark matter annihilation to quarks 
is expected to give a signal in conventional direct detection experiments. This is because the dark matter 
scattering cross section is suppressed by loops of heavy fields, and it is only when dark matter couples directly 
to quarks that the WIMP-nucleon cross section is large enough to give a signal in upcoming experiments. We 
assume in this section that dark matter annihilates predominantly to first generation quarks/lcptons. The 
baryon-number-violating interactions in WIMPy baryogenesis can also induce proton decay due to WIMP 
scattering, but we find that our models arc consistent with all current and projected experimental bounds. 

We first present the Feynman diagrams for the leading processes relevant to direct detection. With WIMP 
annihilation to leptons, we show the diagrams for scattering in direct detection experiments in Figure [121 We 
show the corresponding diagrams with WIMP annihilation to quarks in Figure fT3l 

WIMP annihilation to leptons: X can elastically scatter off electrons at one loop and nucleons at two loops. 
However, direct detection experiments are on the verge of testing dark matter models with nucleon scattering 
at one loop |16[ [51 [5] and electron scattering at tree level (TJ- Therefore, the elastic scattering signals from 
our models are too small to be detected at near- future experiments. 

As we discussed in Section 13. 2[ there can be lepton- number-preserving dark matter annihilation channels 

^Precise reach depends on my,, mg, Ax, A^, and t. 
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in addition to those responsible for baryogenesis, although the parameter space is more restricted in this case. 
The lepton-number-preserving WIMP interactions can lead to conventional WIMP signals in direct detection 
experiments, but do not probe the WIMP's lepton-number-violating couplings, which are the crucial ingredi- 
ents of WIMPy leptogenesis. 

WIMP annihilation to quarks: The dominant contribution to the direct detection cross section is the one-loop 
scattering of dark matter off the right-handed up quark. We estimate the dark mattcr-nuclcon cross section: 



nucleon mass, so rriN. 

We determine the direct detection cross section for the benchmark points in Section 14.21 For the point 
mx = 4.25TeV, = 7.25 TeV, tos = 5 TeV, Ax = 2.7 and As = 4.5, we find ctx-at ~ IxlQ-^" cm^. For the 
point mx = 0.9 TeV, = 1.2 TeV, ms = 1.5 TeV, Ax = 0.22 and As = 2.8, we find (Tx-n « 4 x lO^^^ cm^. 
The current limits from the XENONIOO, CDMS experiments [TSl on dark matter direct detection have a 
minimum bound of ^ 10"''^ cm^ for WIMPs with masses of ~ 50 GeV. The upper limit on the cross section 
for a TeV WIMP is ~ 10^^'^ cm^. We therefore see that the cross sections for our benchmark points are below 
current bounds but are large enough that they can give a signal in upcoming direct detection experiments 
such as XENONIT l^^. We leave a detailed study of the direct detection reach for future work. 

There also exists an inelastic scattering process that converts an up-type quark to two down-type anti- 
quarks, as we show in the right-hand graph in Figure [T51 Such an inelastic process can lead to nucleon decays 
induced by WIMP scattering. The dominant process is X p ^ X n tt"*", along with the corresponding processes 
with strange quark production {b quark production is kinematically suppressed). To avoid conflict with proton 
decay experiments, the induced proton decay rate should satisfy bounds outlined in |21| . Comparing our model 
to the Hylogenesis model in |21j , we find that the operator giving rise to induced proton decay in our model is 
dimension-9 {X'^uddn/ A^), whereas the corresponding Hylogenesis operator is dimension-7. At the hadronic 
level, the Hylogenesis process is 2 — > 2, in contrast with our 2 — > 3 process, which gives our model a relative 
phase space suppression ^ l/(27r)'^. Furthermore, the Hylogenesis model has a dark matter mass ^ ©(GeV), 
while the dark matter mass in WIMPy baryogenesis is mx ^ ©(TeV). As a result, the dark matter number 
density is smaller by a factor of (GeV/TeV) in WIMPy baryogenesis, and the incident flux of dark matter 
particles is suppressed. Taking into account all factors, the induced proton decay rate for WIMPy baryogenesis 
has a ^ (57)^ (§sv) ^ 10~^^ suppression compared to that of Hylogenesis. Since induced proton decay is 
on the verge of current bounds for Hylogenesis models with a heavy scale A ^ TeV, the proton lifetime in 
our model is safely above the current bound, and not within the reach of near- future proton decay experiments. 

To summarize, models with WIMP annihilation to leptons typically predict the absence of a signal in 
conventional dark matter direct detection experiments, while models with WIMP annihilation to quarks have 
WIMP-nucleon cross sections below the current bounds but accessible in upcoming experiments. Baryon/lepton- 
number-preserving WIMP interactions can also give a signal in direct detection experiments, but such models 
have a smaller viable parameter space. WIMP scattering can induce nucleon decay in WIMPy baryogenesis 
models, but the proton decay rate is far lower than current experimental constraints. 

5.1.2 Indirect Detection 

Models of WIMPy baryogenesis have indirect detection prospects similar to those in conventional WIMP 
scenarios because the dark matter relic abundance is symmetric and is established by thermal freeze-out. This 
is in contrast with many asymmetric dark matter models, which typically have suppressed indirect detection 
signals due to the fact that dark matter is largely asymmetric today. The only asymmetric dark matter models 

^Because the masses of all fields running in the loop arc similar in mass, there is no significant mass suppression to the cross 
section from evaluating the loop integral. 




(42) 



where fi is the reduced mass of the dark matter- nucleon systerr|^. In our models, mx ^ miq, where mjv is the 
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with indirect detection signals are those in which the symmetric component is regenerated after WIMP freeze- 
out [5] . In the following summary, we assume that WIMPs annihilate predominantly through the interactions 
that generate the baryon asymmetry. 

We find that indirect detection is most promising with WIMPy baryogenesis with dark matter annihilating 
to quarks. In this scenario, the final states are color-connected quarks and sterile fields rn, and the quarks 
hadronize in the dark matter rest frame. This populates the low-energy anti-deuteron spectrum, leading to 
a clean, low-background signal at GAPS and AMS-02. The mass reach in this scenario is mx ^ 1 TcV. 
Annihilation of dark matter in WIMPy Icptogcnesis also leads to qq production via Higgs decay, but the 
quarks hadronize in the Higgs rest frame. Fewer low-energy antideuterons are produced, and the mass reach 
is only mx S 200 GeV, which is too low for viable models of WIMPy leptogenesis. We give more details below. 

WIMP annihilation to Icptons: The indirect detection signals are energetic neutrinos, positrons, and secondary 
photons from the leptons produced in WIMP annihilations, along with antiprotons and antideuterons (D) from 
ip'^ —i' h + n bb + n. Unfortunately, the dark matter mass of O(TcV) in leptogenesis gives a flux lower than 
the sensitivities of most upcoming indirect direction experiments. With the standard cross section for thermal 
WIMP annihilation (((Tann) ~ 3 x 10^^^ cni-^/s), the reach of most current experiments like Fermi-LAT [52] is 
in the mass range < 0(100 GeV). One exception is in the scenario with a very steep dark matter profile in 
the galactic center, which occurs in halo models favored by hydrodynamical simulations. In this case, HESS 
measurements of gamma rays from the galactic center arc within a factor of two of constraining a 3 TeV WIMP 
with standard annihilation cross section to leptonic final states [5^. Based on the HESS analysis, it is likely 
that with more data Fermi-LAT could rule out WIMPy leptogenesis models with masses < few TcV if the 
similar assumptions on dark matter distribution arc applied. Such constraints sufii'er from large uncertainties 
in the dark matter profile, however, and we caution that such strong limits on WIMP masses may not be 
possible. 

According to the general analysis performed in [24j , the mass reach of low energy antidcutcron detection 
experiments at AMS-02 and GAPS could be up to ~ ITeV if hadronization happens mostly in the rest frame 
of dark matter annihilation, as occurs in the gg channel in |24[ . However, hadronic decay products in the 
leptogenesis scenario are secondary or tertiary, and hadronization typically happens in the boosted frame, 
similar to the WW channel in the same reference. The resulting mass reach could be only ~ 200 GeV, which 
is too low for WIMPy leptogenesis because sphalerons are decoupled during the era of asymmetry generation 
for dark matter masses in this range. 

WIMP annihilation to quarks: The possible signals are p, D, and 7. In contrast with leptogenesis, the baryon 
asymmetry can be generated after the electroweak phase transition and the dark matter mass can be as low as 
^ 290 GeV according to the bound in Section 15.21 This is promising for detection at upcoming experiments, 
particularly low energy anti-deuteron searches. Because the primary products of WIMP annihilation now 
involve color-connected u and ip, a large proportion of hadronization proceeds in the rest frame, resulting in 
a larger rate of D production [53]. This extends the mass reach at GAPS and AMS-02 [55] to ~ 1 TeV and 
covers a large part of the WIMPy baryogenesis parameter space. Higher WIMP mass regions TeV) may be 
constrained by Fermi-LAT gamma ray observations of the galactic center, but as discussed above with WIMPy 
leptogenesis, these constraints are highly dependent on the dark matter profile [23| . 

In general, models of WIMPy baryogenesis and leptogenesis satisfy all current constraints from indirect 
detection experiments, and future searches for antideuterons are promising discovery channels for models with 
WIMP annihilation to quarks. 

5.2 Collider Detection 

We consider the LHC constraints and detection prospects for new charged particles predicted in WIMPy 
baryogenesis. We find that the LHC can strongly constrain the scenario with WIMP annihilation to quarks 
but may not constrain WIMP annihilation to Icptons. Searches for supersymmctry (SUSY) with missing 
energy are relevant to our models, since WIMPy baryogenesis predicts new charged fields decaying to Standard 
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Figure 14: LHC clcctroweak pair production of and its subsequent decay in the model with WIMP annihi- 
lation to leptons. -0° decays to a Higgs boson and the light neutral fermion rif, while ^/j^ decays to and 



Model particles and neutral fermions. We focus on existing LHC searches for SUSY and leave for later work 
the optimization of collider searches for the particular charged fields found in our models. 

The strongest LHC constraints are bounds on new colored fields, such as giuinos and squarks. Current 
searches at the LHC therefore constrain the scenario with WIMP annihilation to quarks, which has new colored 
fields ip. The luminosity at the LHC is not yet large enough to bound electroweak production of new particles, 
due to the smaller production cross section, and softer jets and missing energy. As a result, current collider 
constraints on in the scenario with dark matter annihilation to leptons are well below the range needed 
for viable WIMPy leptogenesis. Higher luminosity and new targeted searches can improve the LHC reach for 
depending on its decay modes. 

We now consider each scenario in more detail. 

WIMP annihilation to leptons: A characteristic feature of the WIMPy leptogenesis model in section [3] is 
the presence of an exotic vectorlike SU(2)l doublet ij). The neutral and charged components of ij) can be 
pair-produced via electroweak gauge bosons. According to our arguments in Section [3. 1.11 if) decays promptly. 

The dominant decay of tp^ is to Higgs -I- through the interactions in (|18p . where is the light neutral 
mass eigenstate after Higgs-induced mixing between ip and n. The resulting collider signature for pair pro- 
duction is 'ip^ij}^ — >• 46(4j) +^T- The charged component, i/"^, decays to + ni, with a collider signature 
for pair production of — >■ W~^W~ + The relevant diagrams are shown in Figure fT4l 

Searches at LEP constrain the masses of the charged and neutral components of tp with bounds on pair 
production of charginos (x^ — >■ W^x^), m^± > 100 GeV [^Hj. decays look identical to chargino decays, 
so TO^± > 100 GeV as well. Hadronic chargino decays, which have a 4j -I- 0t final state, constrain the 
mass. The LEP bound groups hadronic chargino decays with other decay modes, so the bound is not directly 
applicable to ip'^. A more careful analysis (that we leave for future work) is needed to determine the precise 
bound on ■0", but we expect it to be on the order of 100 GeV as well. The bounds on both and tp^ are 
well below the typical required for WIMPy leptogenesis. 

With the current luminosity of 5 fb~^ at \/s = 7 TeV, the LHC bounds the masses of weakly charged 
particles appearing in cascade decays of colored particles, but does not constrain particles such as ip that are 
only produced directly from electroweak gauge bosons. Therefore, the LHC does not bound at present, 
and the LEP constraint remains the most important. Searches for direct chargino and slepton production with 
future LHC data will improve the bounds on to masses on the order of a few hundred GeV, but this is 
still smaller than needed in WIMPy leptogenesis. 

New LHC searches at 14 TeV could possibly yield stronger constraints. For example. If the Higgs mass 
is known, we could require a reconstruction of the Higgs mass among final state jet pairs, greatly reducing 
backgrounds. If ^ rrih, the final state Higgses are boosted and can be studied with jet substructure 
techniques, as suggested in [57]. 

In summary, the collider constraints on ip are currently too weak to place any bounds on WIMPy leptoge- 
nesis models. Future LHC running will improve the bounds on m^, and we have outlined some of the possible 
signals here. A more detailed collider analysis is deferred to later work. 
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Figure 15: LHC pair production of ip and its subsequent decay in the model with WIMP annihilation to 
quarks. 



WIMP annihilation to quarks: As with WIMP annihilation to Icptons, there are new charged states at the 
weak scale. In this model, tp carries color charge, and the bounds are consequently much stronger than 
for WlMPy leptogenesis: > 590 GeV as inferred from the current LHC gluino search at 7 TeV. The 
phenomenology depends on how ip decays, and we outlined two possible models in ([55)1 and ([59]) . In both, ip 
decays to two jets plus a singlet, and so the collider phenomenology is identical. For the purposes of notation 
in this section, we assume that ip decays through an intermediate colored scalar (f> to 2 Standard Model quarks 
and a singlet fermion, n. 

Ip can be pair-produced at the LHC with the signature — > tp-ip — > Aj+J^T- We show the relevant diagram 
in Figure [T51 

Gluino searches at LHC7 bound the tp mass. Both gluinos and tp decay to jj + and the bounds on 
both gluinos and ip are comparable, as their cross sections differ only by a group theory factor. We correct 
the gluino bounds for this factor. We apply simplified model searches from ATLAS, which place bounds on 
gluino and squark masses in the presence of a massless neutralino [28]. The corresponding fields in WlMPy 
baryogcnesis are ip, the colored scalar 0, and the massless singlet fermion n. The lower bound on the ■0 mass is 

> 590 GeV, which comes from the gluino bound when the squark is much heavier than the gluino. In our 
scenario, this means that m0 > 590 GeV when jti^ ^ m0 (numerically, to^ > 1.2 TeV). The bounds on 
cut significantly into the allowed parameter space for dark matter because 2mx > , and so mx ^ 295 GeV 
for heavy (p. The bounds are stronger for lighter (p because (p and ip can be jointly produced. For example, 
the bound on m.^ is about 30% higher for = 1 TcV. 

The LHC search reach for gluinos is expected to be rUg « 1.44 TeV at 100 fb~^ and i/s = 14 TeV [29] 
(with the assumptions of mSUGRA and heavy squarks) , so models of WIMPy baryogcnesis will be strongly 
constrained by future running of the LHC. The LHC will not reach the highest-mass regions of WIMPy 
baryogcnesis, but will exclude models with masses m^, < 2 TeV, mx S 1 TeV, and 0(1) couplings. 

LHC searches also constrain the mass of the colored scalar (p. Since is not directly relevant to the 
outcome of WIMPy baryogcnesis (apart from the requirement that it be light enough for ip decays to be 
in thermal equilibrium), bounds on do not directly constrain WIMPy baryogcnesis. Nevertheless, the 
production rate of (p is comparable to that of squarks and is very high at the LHC. With the interaction 
([38]) . (p decays to di + 0t and has an event topology identical to squark pair production in the MSSM: two 
jets (possibly 6-tagged) plus missing energy. The current model-independent constraint is > 875 GeV for 
degenerate squarks of the first two generations [28]. In WIMPy baryogenesis, however, only a single field (p is 
necessary, so the bound can be relaxed. Since (p can decay into 6, the bound is approximately that of a sbottom 
squark from D0, m^^ > 250 GeV [30]. Future LHC running at 14 TeV will improve the bound to ~ 2 TeV at 
100 fb~^ [29], and has the potential to discover colored scalars in the mass range of WIMPy baryogenesis. 
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Figure 16: A set of two-loop contributions to the electron EDM that vanishes when summed together. 



5.3 Electric Dipole Moment Constraints 

A viable mechanism for baryogenesis necessitates the existence of new CP phases. Bounds on the electron 
and neutron electric dipole moments (EDMs) strongly constrain many new sources of CP-violation, but we 
find CP phases in WIMPy baryogenesis are not constrained by EDM experiments. The minimal models of 
WIMPy baryogenesis presented in this paper couple new fields to either left-handed or right-handed light 
fermions (but not both), resulting in suppressed EDMs that are consistent with current observations. As a 
result, minimal models of WIMPy baryogenesis do not have the CP problem often associated with models of 
weak-scale physics. 

In the models presented in Sections |3] and lU the fields S and tp couple exclusively to cither left-handed or 
right-handed quarks and leptons. As a result, loops contributing to light fermion EDMs are helicity-preserving 
with an even number of Yukawa couplings, half of which are of the form Xai and the other half A^^-. By 
summing over all permutations of different flavors of S, L and ijj on the internal lines, it can be shown that all 
one- and two-loop diagrams appear in pairs that are complex conjugates of one another. Summing over each 
set of pairs leads to a result that is real, and hence a vanishing EDM. As an example, we show in Figure [TBI a 
set of four diagrams contributing to the electron EDM: the sum of the first two is proportional to XLai^^ pi, 
while the sum of the second two is proportional to ^Lai^L i3i- Therefore, the sum of all four is real and does 
not contribute to the electron EDM. 

The two-loop EDM in the Standard Model vanishes for the same reason as in WIMPy baryogenesis: 
CP-violation arises only in couplings to one chirality of fermion. In both the Standard Model and WIMPy 
baryogenesis, the neutron EDM is non-zero at three loops, and we show the relevant diagrams in Figure [T7l 
The principal difference between the two is that CP violation vanishes in the Standard Model with fewer than 
three generations, so the Standard Model EDM is suppressed by mixings involving all three generations. By 
contrast, WIMPy baryogenesis has a contribution to the EDM with only two generations of quarks that couple 
to more than one flavor of S, and if the model is minimally flavor violating, the EDM will be suppressed 
by sin^ 6c ~ 0.05, the square of the Cabibbo angle. The naive estimate for the neutron EDM in WIMPy 
baryogenesis with 0(1) couplings is 

dn sin^ 0c ruu . . 

e ^ (167r2)3 ml' ^ ' 

Substituting ms ^ 5 TeV and mf ^ MeV gives d„/e ^ 5 x lO^'^^ cm, which is well below the current 
experimental bound of d„/e < 2.9 x 10"^^ cm [31]. The electron EDM from WIMPy leptogenesis is even 
smaller than this, as flavor-changing effects in the charged lepton sector are suppressed by neutrino masses, 
and the EDM is also well below the experimental bound of rfg/e < 1.05 x 10^^^ cm [32] . 

Phases from WIMPy leptogenesis can also contribute to EDMs via other new, weak-scale flelds not included 
in our minimal models. Since these contributions arc model-dependent, we do not consider them further. 

Minimal models of WIMPy baryogenesis do not suffer from a CP problem and are consistent with low- 
energy experiments, but it is possible that other, model-dependent contributions to the EDM could be con- 
strained by and give rise to signals in electron and neutron EDM experiments. 
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Figure 17: (left) Three loop EDM in the Standard Model, (right) The analogous diagram for quark EDMs 
in WIMPy baryogenesis. The photon line attaches to any charged internal fields. 



6 Conclusions 

In this paper, we explored a novel scenario called WIMPy baryogenesis that extends the WIMP miracle by 
generating the observed baryon asymmetry through annihilations of weak-scale dark matter and provides a 
dynamical connection between the dark matter and baryon abundances. We found that natural couplings and 
weak-scale masses for the new fields can lead to the correct baryon asymmetry. As a by-product of linking 
baryogenesis with dark matter annihilation, we introduced a new mechanism for weak-scale baryogenesis, 
avoiding any conflicts with reheat bounds in supersymmetrie theories. 

The key observation is that, if dark matter annihilation proceeds via CP— and (Standard Model) B— or 
L— violating operators, all Sakharov conditions for baryogenesis are satisfied. Successful models also suppress 
washout prior to dark matter freeze-out. In our models, such suppression results from the heaviness of the 
field tp carrying Standard Model gauge charges and B~ or number that is one of the final states in dark 
matter annihilation. Additional discrete symmetries forbid such exotic fields from decaying back to Standard 
Model fields. We presented models where dark matter annihilates to either quarks or leptons, and found viable 
parameter spaces with natural couplings and TeV-scale masses in both scenarios. 

In models where dark matter annihilates to leptons, the lepton asymmetry must be generated before the 
electroweak phase transition so that the asymmetry can be transferred to baryons via sphalerons. As a result, 
dark matter and masses must be ©(TeV). Because the new states are heavy and dark matter does not 
couple directly to quarks, dark matter in this set-up is not in reach of near-future direct and indirect detection 
experiments. In this scenario, V' is charged only under weak interactions, making LHC searches challenging, 
although targeted searches at high integrated luminosity may allow discovery. 

If dark matter annihilates to quarks, baryogenesis can occur after the electroweak phase transition, allowing 
smaller dark matter and ■0 masses. With lighter new states and colored objects, dark matter in these models 
can be within reach of direct detection experiments and indirect detection searches for antideuterons. LHC 
searches for are similar to gluino searches and can exclude < 1.44 TeV at 100 fb~^ and \/s = 14 TeV. 

In both scenarios, WIMPy baryogenesis models can generate both the correct dark matter relic density 
and the baryon asymmetry at the weak scale. Such models predict new weak-scale particles that can lead to 
signals in dark matter direct and indirect detection experiments, and that may be accessible at the LHC. 
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